Weak Pion Production off the Nucleon 



E. Hernandez, 1 J. Nieves, 2 and M. Valverde 2 

1 Grupo de Fisica Nuclear, Departamento de Fisica Fundamental e IUFFyM, 
Facultad de Ciencias, E-37008 Salamanca, Spain. 
2 Departamento de Fisica Atomica, Molecular y Nuclear, 
Universidad de Granada, E-18071 Granada, Spain 

We develop a model for the weak pion production off the nucleon, which besides the Delta pole 
mechanism (weak excitation of the A (1232) resonance and its subsequent decay into Nit), includes 
also some background terms required by chiral symmetry. We re-fit the C^(q 2 ) form factor to 
the flux averaged v^p /i^pir + ANL q 2 — differential cross section data, finding a substantially 
smaller contribution of the Delta pole mechanism than traditionally assumed in the literature. 
Within this scheme, we calculate several differential and integrated cross sections, including pion 
angular distributions, induced by neutrinos and antineutrinos and driven both by charged and 
neutral currents. In all cases we find that the background terms produce quite significant effects 
. and that they lead to an overall improved description of the data, as compared to the case where only 

' the Delta pole mechanism is considered. We also show that the interference between the Delta pole 

and the background terms produces parity-violating contributions to the pion angular differential 
cross section, which are intimately linked to T— odd correlations in the contraction between the 
leptonic and hadronic tensors. However, these latter correlations do not imply a genuine violation 
of time reversal invariance because of the existence of strong final state interaction effects. 
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I. INTRODUCTION 



The pion production processes from nucleons and nuclei at intermediate energies are important tools to study 
the hadronic structure and have become very important in the analysis of the neutrino oscillation experiments with 
atmospheric neutrinos. The energy spectrum of atmospheric neutrino at Kamioka jlj is such that the weak pion 
production contributes about 20% of the quasielastic lepton production and it is a major source of uncertainty in 
the identification of electron and muon events. In particular, the Neutral Current (NC) tt° production contributes 
to the background of production while ir^ contributes to the background of fi^ production. This is because both 
particles, i.e. 7r° and e or and ^ produce similar single-ring events in Cherenkov detectors, commonly used 
O ' in neutrino oscillation experiments. Moreover, the neutral current ir° production might play an important role in 
distinguishing between the two oscillation mechanisms — ► v r and — ► f s teriie |2t- These comments apply also for 
• v e appearance experiments such as K2K Q and MiniBooNE 

The neutrino oscillation experiments are generally performed with detectors which use material with nuclei like 
12 C, 16 0, etc. as targets. It is therefore desirable that nuclear medium effects be studied in the production of leptons 
and pions induced by the atmospheric as well as accelerator neutrinos used in these oscillation experiments. To 
this end, the starting point should be a correct understanding of the reaction mechanisms in the free space. In this 
context we study in this work, the weak pion production off the nucleon driven both by Charged Currents (CC) and 
NC at intermediate energies. The model derived here will allow us to extend the results of Refs. 0, Q for CC and 
Ref. Q for NC driven neutrino-nucleus reactions in the quasielastic region, to higher excitation energies above the 
pion production threshold up to the A(1232) peak. 

In the past, there have been several studies of the weak pion production off the nucleon at intermediate energies [1]- 
Most of them describe the pion production process at intermediate energies 1 by means of the Delta pole (AP) 
mechanism (weak excitation of the A(1232) resonance and its subsequent decay into Nn) and do not incorporate any 
background terms. Here, we have also included some background terms, required by chiral symmetry. Starting from a 
SU(2) non-linear a model involving pions and nucleons, which implements the pattern of spontaneous chiral symmetry 
breaking of QCD, we derive the corresponding vector and axial currents (up to order 0(1/ f^)) which determine the 
structure of the chiral non-resonant terms. Some background terms were also considered in Refs. [ill, EH and [h|. 
In the latter reference, the chiral counting was broken to account explicitly for p and lo exchanges in the t— channel, 



Higher resonance effects, which are certainly important for energies larger than those considered in this work, are carefully discussed in 
Ref. 0. 



2 



while the first two works are not consistent with the chiral counting either, since contact terms were not included. 
Moreover in [ll], [l2j a rather small axial mass (« 0.65 GeV) was used. 

We will show that the background terms produce quite significant effects, which will require to re-adjust the C^(q 2 ) 
form-factor that controls the largest term of the A— axial contribution. We will find corrections of the order of 30% 
to the off diagonal Goldberger-Treiman relation when the Argonne bubble chamber cross section [20| are fitted. Such 
corrections would be smaller if the Brookhaven bubble chamber data [2l| were considered. We will also show that 
interference between the AP and the background terms produces parity-violating contributions to the pion angular 
differential cross section, which are intimately linked to T— odd correlations in the contraction between the leptonic 
and hadronic tensors. However, these T— odd correlations do not imply a genuine violation of time reversal invariance 
because of the existence of strong final state interaction effects. 

The paper is organized as follows. After this introduction, in Sect. |TT] the model for CC neutrino and antineutrino 
induced reactions is presented. There, some general definitions involving kinematics and differential cross sections 
are given (Subsect. Ill A[) . The consequences of isospin symmetry are exploited in Subsect. Ill Bt while in the next 
subsection the model for the WN — > N'it reaction is presented. In Sect. the findings of the latter section are 
extended to the case of NC driven processes. Results are presented and discussed in Sect. II VI and the main conclusions 
of this work can be found in Sect. [Vj Besides in Appendix YK[ the cross section dependence on the pion azimuthal 
angle is discussed in terms of Lorentz, parity and time-reversal invariances, and finally in Appendix [B] we discuss in 
some detail the effects on the neutrino and antineutrino induced cross sections of different relative signs between the 
axial and vector W^AN form- factors, and between the resonant and chiral non-resonant contributions. 



II. CC NEUTRINO AND ANTINEUTRINO INDUCED REACTIONS 



A. Kinematics and differential cross section 



We will focus on the neutrino-pion production reaction off the nucleon driven by charged currents, 

n {k) + N(p) r(fe') + N(p') + 7r(k v ) 



(1) 



though the generalization of the obtained expressions to antineutrino induced reactions is straightforward. 

The unpolarized differential cross section, with respect to the outgoing lepton and pion kinematical variables, is 
given in the Laboratory (LAB) frame (kinematics is sketched in Fig[T]) by 2 




FIG. 1: Definition of the different kinematical variables used through this work. 



dn(k')dE'dn(k n ) 



\k'\ G 2 f +x d\K\\I^l {u) M 



(2) 



with k and kl the LAB lepton momenta, E' = (kl 2 + m 2 ) 1 / 2 and mi the energy and the mass of the outgoing lepton 



(m M = 105.65 MeV, 



0.511 MeV ), G = 1.1664 x lO" 11 MeV" 2 , the Fermi constant, k„ and E v = {kl + ml) 1 ' 2 



2 To obtain Eq. (T2]l we have neglected the four-momentum carried out by the intermediate W— boson with respect to its mass (Mw), and 
have used the existing relation between the gauge weak coupling constant, g = e/ sin 8w, and the Fermi constant: G/v2 = g 2 /8M^y, 
with e the electron charge, 0w the Weinberg angle and Mw the W— boson mass. 
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the LAB momentum and energy of the outgoing pion 3 , and L and W the leptonic and hadronic tensors, respectively. 
The leptonic tensor is given by (in our convention, we take £0123 — +1 and the metric = (+, — , — , — )): 

l$ = (4 v) W + i (4 I/) W = fc ^ + fe ^-^ fc - fe/ + ie M^/3 fe/Q ^ ( 3 ) 

and it is not orthogonal to g M even for massless neutrinos, i.e, L^Jq 11 = —m 2 k a . 
The hadronic tensor includes all sort of non-leptonic vertices and it reads 

rtl ^W = — [ d3p ' 1 

spins * 



WcJ M = J J^^^P' + ^-q-p)(N'n\^ c+ (0)\N){N'n\^ c+ (0)\Nr (4) 



with M the nucleon mass 4 , q — k — k' and E' N the energy of the outgoing nucleon. The bar over the sum of initial 
and final spins, denotes the average on the initial ones. As for the one particle states, they are normalized so that 
(pW) — (27r) 3 2po ( 5 3 (p — p')> an d finally for the charged current which couples to the W + we take 

j? c+ = *„7"(1 - 75)(cos 6 c ^d + sin 61c*,) (5) 

with and ^> s quark fields, and 6c the Cabibbo angle (cos^c = 0.974). Note that with all these definitions, 

the matrix element (N' ir\j£ c (0)\N) is dimensionless. After performing the d 3 p' integration, there will still be left an 
energy conserving Dirac's delta function (S(p'° + k® — q° — p )) in the hadronic tensor, which can be used to perform 

the d\hx\ integration in Eq. ([2]). Since the quantity J dQ^L^J (Wgc-k) IS a scaiar to evaluate it we take for convenience 
q in the Z direction. Referring now the pion variables to the outgoing wN pair Center of Mass (CM) frame (as it 
is usual in pion electroproduction) would be readily done by means of a boost in the Z direction. Note that the 
azimuthal angle fa is left unchanged by such a boost. 
By construction, the hadronic tensor accomplishes 

(W&J V) = (W&jW + i (W&J: ] (6) 

with (Wqq )^ and (Wqq^) real symmetric and antisymmetric parts, respectively. 
As it is explicitly shown in Appendix \^ Lorentz invariance restricts the fa dependence, 

- ^P-^o {A + B cos fa +C cos2&r + D sin fa + Esin2fa} (7) 



dn(k')dE'dSl(k„) \k I 47r2 

with A, B, C, D and E real, structure functions, which depend on q 2 , p ■ q, k v ■ q and k^ ■ p. 
For antincutrino induced reactions we have 



:en 

which couple to the W + and W~ bosons, respectively. 



and we will discuss below the existing relation between the matrix elements of j^ c+ and j£ c _ = j^c+j charged currents 



B. Isospin relations 

The non-strange parts of j^ c ± behave as the spherical ±1 component of an isovector, since 



^4 = -%^ g , iM> u = ¥,-^¥ g , %=^j, 7 -i= r „4 1 = z F ^__» (9) 



3 For m„-, we use the isospin averaged pion mass. 

4 We take the average of the neutron and proton masses. 
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with t the Pauli matrices. Thanks to the Wigner-Eckart's theorem, we find that all (N 1 Tr\j^ ci .(Q)\ N) matrix elements 
are just determined by two of them, f.i. {pir + \ j^ c+ (0) \p) and (mr + \j^ c+ (0)\n) , 

( P n°W c+ (0)\n) = ~[<pK + \j^(0)]p)-(nir + \j^(p)\n)] (10) 

(P7r-|^ c _(0)b) = (nn+\j? c+ (0)\n) (11) 
(nw-\j^_(0)\n) = (p7r+|j£ + (0)|p) (12) 

(n7r°|j£_(0)|p) = -(pn°\j» c+ (0)\n) = -L [(pw+\j^ c+ (0)\p) - (n7r+|j£ + (0)|n)] (13) 

Thus, Eqs. ([5]) and (fTT]) " (p~5|) allow us to determine all CC antineutrino cross sections from the neutrino induced 
amplitudes. 

Besides, the vector contribution of the matrix elements of the weak CC between N and N'it states is related to 
that of the electromagnetic current s£,' m (0), 

= \^ q l^ q -\^sl^s + ^f q l^^ q (15) 

The matrix elements of the isovector part (tq) are related to those of j C c±, while the first two terms are isoscalar 
operators. One easily finds, 

-J— (pTr+l^^b) = V2{nn°\sU0)\n) + (pK-K m (0)\n) (16) 

COS VQ 

-J—(nn+\V» + (0)\n) = ^(p^KmWb) - (prr"^^) (17) 

COS (7c 

C. Model for the WiV -» iV'vr reaction 
i. SU(2) non-linear a model 

Let us start with the effective lagrangian of the SU(2) non-linear a model. It implements the pattern of spontaneous 
chiral symmetry breaking of QCD and it is given by 

C Nv = %vf [0„ + V M ] * - M§* + 5A * 7 ^75^* + ^Tr [c^cW] (18) 
where ^> = ( ^ ] is the nucleon field. The fields V M and are given in terms of the matrix field £ derived from the 



n 



pion fields 5 



V„ = g + ^ = i - ^) (19) 



The pions (j> are the Goldstone bosons associated to the spontaneous breaking of the SU(2)yxSU(2)yi chiral sym- 
metry We describe their dynamics in terms of 2 x 2 matrix field U given by 

V 2 v 2 

with ~ 93 MeV the pion weak decay constant. The matrix field 

£ = e W/(2/„) (21) 



5 We use a convention such that <j> = {4> x — i<j>y)/V2 creates a n from the vacuum or annihilates a n + and the <f> z field creates or 
annihilates a 7r°. 
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transforms under SU(2)yX 811(2)^4 as 

SU(2) y 



SXJ(2) A 

T v l;Tl £ =J- TkAt=A^T 



(22) 



where Ty = exp(— if • By)j2 and T4 = exp(— if • 9a)/2 are global transformations belonging to SU(2)y and SU(2)^ 
respectively. As for A = exp(— if • 9\)/2, it is a unitary matrix field that depends on the axial transformation Ta and 
the <j> Goldstone boson fields. 

On the other hand the nucleon field, \& transforms as 



SU(2) y 



SXJ(2) A 



(23) 



Each term of the effective Lagrangian of Eq. (fT8|) is separately invariant under the chiral group SU(2)y xSU(2),4. 
This is why one can introduce an axial nucleon coupling gA I in the model without violating chiral symmetry. We 
will use qa — 1-26 through this work. 

Explicit SU(2)^ breaking terms are included in the model as 



ml^Tr(U + tf-V2f n ) 



to give mass to the pions. Neglecting 0(1/ f^), the effective Lagrangian of Eqs. (fT5|) and (|24p reads 



C = *[i^-M]* + i9 M 0^0-im^ 2 



ra _ 9a 

Mnt f 

J; 



— t 7" — t — t — * 7~ — > 

(f) 2 -d^- {<t>d^)-(j) 



(24) 



(25) 



(26) 



In contrast to the linear cr model, the coupling between the p and n and the pions is of the pseudo-vector type. 
Writing the coupling constant in the usual way as g 7t Nff/2M — f/m n , we recover the Goldberger-Treiman relation 



, run 
1 = 2TJ A 



(27) 



which phenomenologically is satisfied at the level of 5%. 

The vector and axial currents that we derive from the lagrangian in Eq. (|f 8p and the transformation properties of 



6 



the fields 6 are given by 



. 2 2j7r 



3/5 



* &>4>) +o(- 

Jtt 



(31) 



A' 



(32) 



and determine the weak transition vertex where the W— boson is absorbed. This is because these currents, up to a 
factor of cos 9c, are the hadronic realization of the electroweak quark current jg c for a system of interacting pions and 
nucleons. Thus, A£ and V<j* account for the W— decay into one and two pions, respectively, while and provide 
the WNN vector and axial vector couplings. Besides, A% and V£ lead to contact WNNir vertices and finally 

and V£ either contribute to processes with more than one pion in the final state or provide loop corrections to the 
leading order amplitude for one pion production. 

The overall normalization is fixed by the W + np vertex 



(j>;p'=p + q |j c a c+ (0)|n;j?) =cos0 c u(j>')(V£(q) - A%{q))u(p) 

where the u's are Dirac spinors for the neutron and proton, normalized such that uu 
nucleon currents are given by 



(33) 

2M, and vector and axial 



F 2 v (g 2 
2M 



-a q v 



A a N (q) = G A (q 2 ) x 7 Q 75 + 



79 75 



(34) 



being all form factors real thanks to invariance under time reversal. Invariance under G-parity has been assumed to 
discard a term of the form (p^ + p ,fJ, )"f5 and we have only considered the pion pole-contribution in the pseudoscalar 
form factor. Isospin symmetry relates the vector form factors to the electromagnetic ones 7 



W) = £ Wfl 9 ) - *!V)) 



and the axial form-factor is given by [23j 

G A (q 2 ) 



9A 



[1-ql/Mir 



HnF?(q 2 )) 



g A = 1.26, M A = 1.05 GeV 



(36) 



(37) 



For infinitesimal vector and axial transformations, we get from Eq. J22D 



■e A ='^ r L e A + o(^) (28) 

Jit 



SU(2) V . SU(2) A ^ 

&$ = (0 V x $) + o(i//|) &$ = ue A + (^{<i>-o A )-8 A <t>' 2 ) /(3/x) + o(i//|) 

SU(2) V SU(2) A 

where 8<j> and denote the infinitesimal variations of the fields. 
7 We use the parameterization of Galster and collaborators (2 2f | 



Up Tm, f „ Ar U n 7 "Up U n7 U n 7 . s Up 

'1 " 



V ' L y UnF? = ™ E , G V E = = = -(1 + X n r)^L = ( l—y. ) (35) 



1 + T ^ 1 + T fi fi p fJ, n VnT \l-q 2 /M 2 D 

with r = -g 2 /4M 2 , M D = 0.843 GeV, n P = 2.792847, /j, n = -1.913043 and A n = 5.6 



(29) 
(30) 
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Thus, one realizes that, up to an overall — \/2cos8c factor, the +1 spherical component 8 (L4^]+i = 
— (L4^] x + iL4^] y ) ) of A% gives the axial vector contribution, at q 2 = 0, of the W + n — > p weak transition. Be- 
sides, — \/2[^4a]+i and the ttNN coupling in Cf nt lead to the q Mr ys term in Eq. (l3"4")) . Analogously, — V2 cos 9c[V^ ]+i, 
provides the Dirac part of the vector contribution, at q 2 = 0, of the W + n — > p weak transition. The magnetic part in 
Eq. (|34p is not provided by the non-linear sigma model constructed here which assumes structureless nucleons. 

From the above discussion, we conclude that — v^cos^c ([V^+i — provides the W + — absorption vertex, 

with the appropriate normalization, in the (N'w\j^ c+ {0)\N) matrix element. We will improve on that by including 
the q 2 dependence induced by the form factors in Eq. and adding the magnetic contribution (F 2 y term) to the 
vector part of the W + N — > N amplitude. 



2. The WNA and NAtt vertices 



At intermediate energies, the weak excitation of the A(1232) resonance and its subsequent decay into Ntt domi- 
nates the WN — » N'tt reaction. A convenient parameterization for the W + n — > A + hadron matrix element is the 
following Q 



(A+; PA = p + q |j£+(0)|n;p) = u a (p A )T a » (p, q) u(p) cos 9 C , 



where 



(39) 



r^(p,g) 



(9^4 - q a r) + jfe (g a »q ■ pa - <z> A ) + ^ (g^q ■ v - gV) + C%<T v. 

% {9 a H - q a r) + % (g a "q ■ PA - q a p A ) + C£g"» + §£gV| , Pa = P + q 



(40) 



with C^4 5 6 an d C3 4 5 6 scalar and real vector and axial form factors, which depend on q 2 . Besides, u a is a Rarita 
Schwinger spinor for the A + . The A^A7T coupling is given by 



C.N A = ^l^ft(^0)* + h. 



(41) 



where ^ u is a Rarita Schwinger .r = 3/2+ field, ft is the isospin transition operator 9 from isospin 1/2 to 3/2 and 
/* = 2.13 x / = 2.14. The Goldberger-Treiman relation implies here 10 



2 U 



c 5 A (o) = Vi— r ~i.i5 

V 3 

For the A— propagator G fll/ (p A ), we use in momentum space 

P^(pa) 



G^(pa) 



p\ - Ml + iM A T A 

with M A the resonance mass (~ 1232 MeV), P^ the spin 3/2 projection operator 



P^{pa) = -0>a + M a ) 
and r a the resonance width in its rest frame, given by 



aT - - - 



zp A p A , ip A r-Pir 



3 Ml 3 



Ma 



(42) 



(43) 



(44) 



67T Vm.r/ V s 



A2(s,m2,M 2 ) 



e(Vs-M 



s =Pa 



(45) 



Note, 

- V2r+i\n) = (t x +ir y )\n) =2\p) (38) 

It is a vector under isospin rotations and its Wigner-Eckart irreducible matrix clement is taken to be one. 
Note that the sign is quoted incorrectly in Ref. (see comment in Ref. |ld|). 
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with X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz and 9 the step function, as deduced from the lagrangian of Eq. (|4ip . 

The determination of the form factors follows from general principles and experimental results. The imposition of 
the conserved vector current hypothesis implies C$ — 0. The other three vector form factors are then given in terms of 
the isovector electromagnetic form factors in thep — A + transition. The analysis of photo and electroproduction data 
of A is done in terms of multipole amplitudes [2J] . Most of the previous papers on neutrino production [l3|, [l?], HH 
assume Mx+ dominance 11 of the electroproduction amplitude, which implies C\ — and a relation between C\ and 
C3 . Here, we take advantage of the recent work of Lalakulich, Paschos and Piranishvili and improve on that by 
including the effect of the subdominant multipoles [Tg} 



2.13 



(1 



f/M 2 ] 



4M? 



cl = 



-1.51 



(1 



q 2 /M 2 ) 2 



q- 



0.48 



(1 



f/M 2 ) 2 



0.776M?, 



(46) 



with My = 0.84 GeV. Among the axial form factors the most important contribution comes from C§ whose numerical 
value is related to the pseudoscalar form factor C§ by PCAC. Since there are no other theoretical constraints for 
C^C? 2 )? C^(q 2 ) and C§(q 2 )/C^(0), they have to be fitted to neutrino scattering data. The available information 
comes mainly from two bubble chamber experiments, ANL [2(| HH, H3 an d BNL [2l|, [28[. We adopt the Adler model 
H, as the ANL and BNL analyses did, where 



Ci(q 2 ) = - 



Cf(q 2 ) = 



(47) 



and for C^ 6 we shall use [Tf 



1.2 



1 2 /M 2 AA ) 2 



1 - 



1 — 

A A 



C£{q 2 ) = C£{q 2 ) 



M 2 



withM^A = 1.05 GeV. 



(48) 



The value for was found to give a small contribution to the cross section and set ting the form factor to zero 
is consistent with early dispersion calculations [1, [2i| and recent lattice QCD results [30j . Note that the contribution 
to the differential cross section of the Cq term will be always proportional to the outgoing lepton mass. 

Isospin symmetry implies that the transition matrix element W + p — > A ++ is a factor y3 bigger than the W + n — » 
A + one discussed above. 

For the weak transition A — > N , we have 



(N;p'\^ c+ (0)\A;p A =p'-q) = (A;p A = p' - q\f cc _{Q)\N,P') 

_ (s> 1) §l*iv, -1>*a) 



(I 1 3| _ 1 1 I) 

\ 2 ' ' 2 I 2 ' ' 2 / 



{u a (p A =p'~ q)T a » (p', -q) u(p' ) cosM 



(49) 
(50) 



with (tx, 1%, t\m\, m^, m) Clebsch-Gordan coefficients and tjv and 4 a t ne nucleon and delta isobar isospin third com- 
ponents, respectively. 



3. Explicit expressions for the {pn + \ (0) |p) and (mr + \j^ c+ (0)\n) amplitudes 

In this subsection, we give explicit expressions for the (pTr + \j^ c+ (0)\p) and (mr + \j^ c+ (0)\n) amplitudes, which we 
will denote by {jc C +)pit+ an d (jc c+)m r+ , respectively. All (7V'7r|j^ c± (0)| N) matrix elements can be expressed in terms of 
these two amplitudes (Subsect. Ill B[) . The model consists of seven Feynman diagrams, depicted in Fig. O constructed 
out of the W+N -> N, W+N -> A, W+N -> iVvr and the contact W+tt -> tt weak transition vertices (E qs. (|3TH33|) 
and p9]) ) and the nNN, ttttN N (Eq. l[2"6"j)) and nNA (Eq. ijil])) couplings, discussed in the Subsects. HTCTI and 
III C 21 Since we have included a q 2 dependence (F^ (q 2 )) on the Dirac part of the vector WNN vertex and to preserve 
vector current conservation, we also include form-factors (Fpp(q 2 ) and FQ T (q 2 )) in the V£ and weak operators. 
This partially accounts for the nucleon structure. 

We do not include loop corrections. This model is an extension of that developed in Ref. [3l[ for the eN — > e'Nir 
reaction, though there exist some minor differences related to the used form factors and a non-relativistic reduction 
was performed in Ref. [3l| . 



11 Recent data determine the contribution from the electric multipole E±+ /M x + to be ~ —2.5% and from the scalar multipole S x + /M±+ ~ 
-2.5% S3 . 
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N 



FIG. 2: Model for the W + N — > N'ir reaction. It consists of seven diagrams: Direct and crossed A(1232)— (first row) and nucleon (second 
row) pole terms, contact and pion pole contribution (third row) and finally the pion-in-flight term. Throughout this work, we will denote 
these contributions by: AP, CAP, NP, CNP, CT, PP and PF, respectively. The circle in the diagrams stands for the weak transition 
vertex. 



The amplitudes read, 



AP 



iC* A — v^cosi 



<c 



p\ - M% + iM A T A 



■u(p')P a p(pA)r 0ll (p, q)u(p), PA=p + q, 



and C*=( *» Pn+ + 
1/3 tor nir^ 



Jcc+ 



;m<1 cCA =(l for Z+)> r^( P ', 9 ) = 7 °[r^(/ ! -g)]V 



CAP 



Jcc+ 



NP 



ACNP 7k C ° S6C U{r)hl \p + q V-M- + l e m<!) AM] ° 



■NP 



for pix^ 

1 for nix" 



Jcc+ 



CNP 



Jcc+ 



9 M 



CT 



PP 



PF 



= -iC 



cl- 



ip' ~ q) 2 ~ M 2 + ie 
-jL- cos 9c u{p')r (9AF^ T (q 2 h 5 ~ F p ((q - fc.) 2 )) u(p), C CT = ( _\ 



CT ( 1 f° r P 71 "" 1 

for 717T" 



-\C PP F P ((q - kn f)^-co S 9c -j^u(p') i u(p), C pp = 1 



V2U 



pp ( 1 for p7r + 
for nir + 



■ lC PF F PF (q 2 )^rCo S 0c J 2K -J r ^ 2Mu(p') l5 u(P), C PF = ( _J 



pp ^ 1 for 

for mr~ 



(51) 



10 



Note that in the PF (PP) term the weak transition is purely driven by the vector (axial) CC. The contribution 
proportional to gA in the CT diagram is the one due to the vector weak transition. We impose 

F PF (q 2 ) = Fct(q 2 ) = 2iW) = F» ~ FT (52) 

to preserve conservation of vector current, as we required for the AN weak transition. Besides, we have included a 
form factor 

Fpit) = ^-r-z, m p = 0.7758 GeV (53) 

1 — t/rrip 

in the PP term to account for the p— meson dominance of the ttttNN coupling. To preserve PCAC, the same form- 
factor has been included in the CT axial contribution. 

In the pion in flight term, we have the coupling ttNN with a virtual pion. It is usual in the literature to use 
a form factor to account for the off-shellness of the pion. To preserve vector current conservation, if one includes 
this form-factor in this term, one should also multiply by the same factor the CT term and the NP and CNP 
contributions. This was the adopted scheme in the study of the eN — > e'Nir reaction carried out in Ref. [3l|, where 
the induced changes by its inclusion turned out to be moderately small. In the weak pion production case, there are 
more important sources of uncertainties 12 and the existing measurements are poorer. Being this work one of the first 
studies of weak pion production where background terms are added to the dominant A contribution, and for the sake 
of simplicity, we do not include this form factor. 

The average and sum over the initial and final spins in Eq. Q is readily done thanks to 

^u{p')S^u{p) [u{p')S°u{p)\* = ^Tr((/ +M)S»U>+ Af) 7 5 t " 7 ) (54) 

spins 

where the spin dependence of the Dirac's spinors is understood and 6> M is a matrix in the Dirac's space for each value 
of the Lorentz index ji. 

III. NC NEUTRINO AND ANTINEUTRINO INDUCED REACTIONS 

The unpolarized differential cross section in the LAB frame for the process 

n {k) + N(p) n{k') + N(p') + (55) 

reads 



|/; '' l ^ dk ^m{W^p (56) 



dO.{k')dE'dn{K) \k | !6tt 2 Jo E„ 

with k 1 , Fj' = \k' | the LAB outgoing neutrino momentum and energy. The leptonic tensor is given in Eq. ^ and it 

is now orthogonal to = (k — k')^ for massless neutrinos, i.e, L^Jq^ = L^q" = 0. 
The hadronic tensor reads 

«cJ M = 7^£ / ^^\p' + K~ q -p)(N'^M\m(N^\f n M\NY (57) 

spins N 

where the neutral current at the quark level is 

j£ c = * u7 ^(l - | sin 2 W - 7 5 )^ u - ^7^(1 - - sin 2 9 W - j 5 )^ d - ^7^(1 - - sin 2 W - 75 )^ s 

- * g 7' I (l-75)r 1 * g -4sin 2 ^^ m -* s 7^(l-7 5 )* s (58) 



12 For instance: partial knowledge of the A resonance form factors, possible pion off-shell effects in the weak transition vertex of the PP 
diagram, etc. . . 
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where 6\y is the Weinberg angle (sin 2 6yy — 0.231). 

Both, lepton and hadron tensors are independent of the neutrino lepton family, and therefore the cross section for 
the reaction of Eq. (|55|) is the same for electron, muon or tau incident neutrinos. For antineutrino induced reactions 
we have, besides the relation of Eq. ([5]) for the leptonic tensor, 

WcJ (P) = (59) 

As discussed above for CC induced process, Lorentz invariance here also restricts the </>„■ dependence, and the NC 
differential cross section can be written as in Eq. ([7]). This (j)^ dependence has been carefully studied in Ref. [HJ as a 
potential tool to distinguish r— neutrinos from antineutrinos, below the r— production threshold, but above the pion 
production one. 

The NC can be expressed as 

j£ = * 9 7 M (l-2sin 2 ^- 75 )T 1 * g -4sin 2 ^< n7S -^7 p (l-75)* s (60) 
where the isoscalar part of the electromagnetic current is given by 

*L,is = \%-f% - l^sl^s (61) 

Isospin symmetry relates the matrix elements of the isovector part (tq term) of jg c with those of the vector (V^+) 
and axial (^4c C+ ) part of the current j£ c+ (= V^ c+ — A^ c+ ), 



(p7r°|* g7 ^(l - 2 sin 2 % - l5 )T^ q \p) = 1 {(1 - 2 sin 2 B w ) \V£ + {0)\p) + (nir + \V£ + (0)\n)] 

v2cost'c* 1 

- [(pn+\A» cc+ (0)\p) + (nn+\A» cc+ (0)\n)] } (62) 
(nn+\* q ^(l -2sm 2 8 w - j 5 )r^ q \p) = L-{(1 - 2sm 2 6 w ) [(pir+ \V» + (0)\p) - (nw + \V£ + (0)\n)] 

COS (7(7 I. 

- [(pn + \A» cc+ (0)\p) - (nn+\A» cc+ (0)\n)} } (63) 
{rm°\%'f{l - 2sin 2 9 W - J 5 )r^ q \n) = ( P 7r°|§ 9 7^(l - 2sin 2 6 W - 75)^*, |p) 

(p 7 r-|* g 7^(l-2sin 2 ^-7 5 )T 1 * g |n) = -(mr + \^ q ^(l - 2 sin 2 W ~ l5 )r^ q \p) (64) 
For the isoscalar part of the electromagnetic current we have 

( n7T+ \ S cmJs\p) = (P^WLj S \ n ) =^(P n °\ S en l Js\p) ^-^( n7r °\ S ernJs\ n ) ( fi 5) 



with 

01 U I \ _ ■"..,;,"".»/ -/'" .--..in 



(p*° I °L is\p) = - ( Tt7r °l s e m (o)l")-^ l^ m (o)b) (66) 



To compute (Nw°\s% m (0)\N) , we derive the electromagnetic current associated to the lagrangian of the SU(2) non- 
linear a model of Eq. (|25|). 

s£ m = ¥7* (^Y^j * + ^7^75 (t\^ + rlrf) * + i (tfd^t - 4>d^ ] ) +■■■ (67) 

where we have only kept those terms contributing to one pion production in absence of chiral loop corrections. Thus, 
within our framework, the model for the 7-/V — > ttN reaction would consist of direct and crossed nucleon pole, contact 
and pion-in-flight terms. As we did for the CC driven processes, such model should be supplemented by including i) 
the q 2 dependence induced by the Dirac Ff' n form factors, ii) the magnetic contribution in the ^NN vertex and iii) 
the direct and crossed A(1232) pole terms [31(. However, the A— resonance diagrams cannot contribute to the matrix 
elements of the isoscalar part of the electromagnetic current. Besides, from Eq. (f6"7|) we see that neither the pion-in- 
flight nor the contact terms contribute for ir° photoproduction. Hence, to compute (rnr°\s£ nl (0)\n) — (p7r°|s^ m (0)|p), 
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we are just left with the direct and crossed nucleoli pole terms 



(0)|») - (p7r°|<C(0)|p) 



M 



(p + g) 2 - M 2 + ie 



/ - + M 
(p' -q) 2 - M 2 




with 



F( s (q 2 ) = § (if (g 2 



•JTGf 8 )) 



MsFl S {q 



(n p F${q*)+H n F?(q 2 )) 



(68) 



(69) 



Finally, we pay attention to the matrix elements of the isoscalar operator ^^7^(1 — 75)^5 which are sensitive to 
the strange content of the hadrons. Due to its isoscalar character we have 



(nvr+l - 75)^) (0)|p) = <p7r-| - 75 )* s ) (0)|n) = V2( p7 r°| (^(1 - 75 )tf s ) (0)|p) 

= -V2(nvr | - 7b)*.) (0)|n) 



(70) 



This part of the NC operator can neither lead to iVA transitions nor couple to a single pion. Thus, and assuming 
a model for the Z°N — > N'tt reactions similar to that used for the CC driven process, we should consider the 
contributions of the direct and crossed nucleon pole, the contact and the pion-in-flight terms to the matrix clement 
of the , I' S 7 M (1 — 75) Is quark operator. The contribution of the first two terms (NP and CNP) reads 



(jm°\ (* s 7"(l-7 5 )* s )(0)|p) 



"i ^75 



i> + i + M 



(p + q) 



-G s A (q 2 )n 5 ~ G s p q^5 
-G s A (q 2 )n 5 - Gfxf 75' 




(71) 



where Ff,fj, s F2,G A and G S P are the strange vector and axial nucleon form factors [33(. The pseudoscalar part of the 
axial current does not contribute to the differential cross section for massless neutrinos and for the rest of strange 
form factors we use the results of the fit II of Ref. [3J| , 



G A {q 2 ) 



as 



{l-q 2 /{M A ) 2 f 



Fl{q 2 )=^{q 2 







(72) 



with g s = -0.15 and M% = M A . 

The vector part in Eq. (|7"Tj) is conserved, i.e., it is orthogonal to g M independently of F{. Because of parity and 
angular momentum conservation, a pion-in-flight type term can only contribute to the vector part of the matrix 
element of the , I , S 7 M (1 — 7s) v I / s operator and its contribution should be proportional to (2k n — q)^u(p')^u(p), as in 
Eq (|51[) . Assuming a structure of the type w(p')757 M u(p ), as in Eq (|51[) . for the contact term vector contribution 
to the matrix element of the ^ s7 M (l — 75)^3 operator, we will conclude that both types of contributions should be 
exactly zero to preserve vector current conservation. Within our scheme, we cannot say anything about a possible 
contact term axial contribution to (p7r°| (\& ^7^(1 — 7s)^ s ) (0) |p) , that for simplicity we will neglect. Thus, we will 
assume that this latter matrix element is given by the NP and CNP contributions in Eq. (fTTj) . 



IV. RESULTS 



In this section, we will show differential and partially integrated neutrino and antineutrino cross sections for pion 
production processes driven by both CC and NC. As it is usual in pion electroproduction, we will work with angular 
pion variables (dn*(k n )) defined in the outgoing nN pair CM frame, while the incoming and outgoing lepton variables 
will be in the LAB frame. Wc will pay here an special attention to the CC pion production cross section dependence 
on the azimuthal pion angle 0* (note that this angle is not affected by the LAB^CM boost), and thus we will show 
the different contributions to S* (q 2 ,p- q, 6* ,</>*), defined by its relation to the differential cross section, 
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d 5 <Tv,i = (73) 

dn(k')dE'dQ*(k^) \k | 4tt 2 

£* (q 2 ,p ■ q,0*, </>*) = {A* + B* cos 0* + C* cos 20* + 13* sin 0* + 12* sin 20* } (74) 

For the NC case this dependence has been already discussed in Ref. [13] with the aim of distinguishing between v r 
and v T below the t— production threshold, but above the pion production one. 

After integrating the pion solid angle, we will take as independent variables the incoming neutrino energy E = \k\, 
the invariant mass W of the outgoing pion-nucleon pair (W 2 = (p + q) 2 ) and the squared of the lepton four-momentum 
transfer q 2 , 



ttW 



dq 2 dW <m(k')dE' ME\k' 



a = i/il, vv. (75) 



where W varies in the range (m^ + M) < W < (V~S — mi), with S = (k + p) 2 = M{M + 2E). Thus, the incoming 
neutrino energy in the LAB system, E, should be greater than (m w + mi + (m, + mi) 2 /2M) for the pion production 
process to take place. Besides for a given outgoing irN invariant mass W, q 2 is comprised in the interval 



,{W) 



2E, 



CM 



E' 



CM 



E' 2 

CM 



< q < 



2E, 



CM 



E' 



CM 



E" 2 

^CM 



(76) 

with Ecm = (S — M 2 )/2\> r S and E' CM = (S — W 2 + m 2 )/2\/~S, the incoming neutrino and outgoing lepton energies 
in the neutrino-nucleon CM frame. It is also useful to perform the phase space integrals in the other order around, 
which allows one to find the da/dq 2 differential cross section. The total range of q 2 is given by q 2 ^ in (W = m^ + M) < 



< 



,(W = jtItt + M) and for a given q , the outgoing ttN invariant mass W varies 



Wn 



M + mn <W < 



M 



2M 



;E- 



AE 



= W max (g 2 ) 



(77) 



where in all equations above, mi should be set to zero for NC processes. 
In the outgoing ttN CM frame we will also use that 



■6(p' + k°-q -p°) = 



\kJE', 



N 



w 



CM 



W 2 + M 2 - mj 



x \i(W 2 ,M 2 ,ml). 



(78) 



Since the main dynamical ingredients of our model are the excitation of the A resonance and the non-resonant 
contributions deduced from the leading SU(2) non-linear a lagrangian involving pions and nucleons, we will concentrate 
in the M + m n < W < 1.3 — 1.4 GeV region. For larger invariant masses, the chiral expansion will not work, or at 
least the lowest order used here will not be sufficient (3f| [3(| . Moreover, the effect of heavier resonances will become 
much more important [19]. Thus, we will limit the available phase-space to guarantee that the invariant mass W 
will lie in the above range. For a fixed incoming neutrino energy, imposing an upper limit in W will lead to different 
amounts of phase-space reduction depending on q 2 (see Fig. [3]). For neutrino energies of about 1 GeV of relevance in 
the CC ANL [2(| and BNL [2l|, [28| bubble chamber experiments, great part of the available phase space satisfies the 
W < 1.3 — 1.4 constrain. As the neutrino energy increases, the q 2 interval which leads to CM irN energies around the 
A— resonance pole gets reduced and the corresponding kinematic cuts performed by the various experiments produces 
a significant reduction of statistics. 

We will see, as it also happens in the pion electroproduction case |3lj . that the inclusion of non-resonant terms 
(background terms) plays a crucial role close to the M + m n threshold. 



A. CC pion production cross sections 



There exist several sets of data taken and analyzed in the late seventies and early eighties. The most detailed 
studies, including measurements, not only of the totally integrated neutrino cross sections, but also of the neutrino 
flux-averaged q 2 and some angular distributions were made in the ANL 12-foot bubble chamber [2(| HtJ and in the 
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FIG. 3: Upper integration limit lV ma x as a function of q 2 (Eq. i77l ) for incoming neutrino LAB energies E = 0.5, 1.0, 1.5 and 5 GeV and 
different outgoing lepton masses. In all case the horizontal line stands for the phase space threshold M + m^. 



BNL 7- foot deuterium- filled bubble chamber [U [28[ . In both experiments the bubble chambers were exposed to a 
wide-band of muon-type- neutrino beams with average energies of approximately 1 GeV (ANL) and 1.6 GeV (BNL) 
and events for the v^p — > [i~pir + , v^n — > [j,~pir and v^n — > n~nir + reactions, with and without a W < 1.4 GeV 
cut, were obtained. The ANL experiment used hydrogen and deuterium targets, though most of data come from 
deuterium exposure. Incoming neutrino energy distributions can be found in Figure 8 of Ref. [37| and in Figure 7 of 
Ref. 38] for the ANL and BNL experiments, respectively 13 . 

Muon-type-antineutrino (energy beam peaked around 1.5 GeV) induced total cross sections off the proton (neutron) 
for final ir~p and ir°n (ir~n) channels with and without the invariant mass cut W < 1.4 GeV were measured in the 
Gargamelle propane experiment at CERN PS [3^ |. 

There also exist experiments at higher neutrino energies carried out at the FNAL 15 foot bubble chamber 
(neutrino energies from 5 to 100 GeV) and at CERN [4l| . In this latter case, a hydrogen target was illuminated with 
a wide band neutrino and antineutrino beams (energies from 5 to 120 GeV), the mean event energy being about 25 
GeV. At such high energies the integrated cross section remains constant with high accuracy, so the exact value of 
neutrino energy is not important. The implementation of the invariant mass cut W < 1.4 GeV reduces significantly 
the statistics and we will not consider these data sets in this work. 

We start looking at the flux averaged q 2 differential cross sections for the reaction v^p — ► [i~pn + measured by 
the ANL and BNL experiments (Fig. In the latter experiment, the cross section overall normalization is not 



Flux— averaged q 2 differential cross sections read 

and similarly for other flux-averaged differential cross sections. 
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FIG. 4: Flux averaged q 2 -differential v^p -> /i~pir+ cross section *° oV dW compared with the ANL [H (left) and BNL [2l[ 

(right) experiments. Dashed lines stand for the contribution of the excitation of the A++ resonance and its subsequent decay (AP 
mechanism) with (0) = 1.2 and M A & = 1.05 GeV. Dashed-dotted and central solid lines are obtained when the full model of Fig. [2] 
is considered with C^(0) = 1.2, M A A = 1-05 GeV (dashed-dotted) and with our best fit parameters C^(0) = 0.867, M AA = 0.985 GeV 
(solid). In addition, we also show the 68% CL bands (solid lines) deduced from the Gaussian correlated errors quoted in Eq. {50}. 



provided. The A vector form factors are fixed by electroproduction data, while the axial weak AN transition form 
factors have been adjusted in such a way that the A ++ contribution alone would lead to a reasonable description 
of the shape of the BNL data (see for instance Ref. [IB]). Moreover, this set of axial form factors also leads to a 
reasonable description [l7| of the ANL data (dashed line in left panel of Fig. 2]). The agreement with the ANL data 
is certainly worsened when the background terms, required by chiral symmetry, are considered (dashed-dotted line). 
Since the C§{q 2 ) controls the largest term of the axial contribution, this strongly suggests the re-adjustment of this 
form-factor. Assuming the same q dependence as in Eq. US), a x 2 — fit to the flux averaged (W < 1.4 GeV) ANL 
v^p — » [i~pir + q 2 differential cross section provides 

C*^(0) = 0.867 ± 0.075, M AA = 0.985 ± 0.082 GeV, (80) 

with a Gaussian correlation coefficient r = —0.85 and a x 2 /dof = 0.4. This fitted axial mass in the weak iVA vertex 
is in good agreement with the estimates of about 0.95 GeV and 0.84 GeV given in the original ANL reference [2(| 
and in the work of Ref. [TtJ • On the other hand, we observe a correction of the order of 30% to the off diagonal 
Goldberger-Treiman relation (Eq. (02])). The lattice QCD results shown in Figure 4 of Ref. [30( might support the 

ratio y/f^.T/C^O) beco ming significantly larger than unity for realistic pion masses. 

The ANL data come mostly from deuterium exposure, and thus deuteron structure effects might affect/are included 



in this determination of the C5 form-factor. Such effects were investigated in Ref. [14j, where it was estimated that 



they were always less than 7%. The solid line in the left panel of Fig. [4] shows the quality of the fit. We also plot 
the 68% CL band deduced from the statistical errors quoted in Eq. (|80|) . We do not fit to the BNL q 2 differential 
cross section because this data set is given in arbitrary units. To fix the overall BNL data scale, we normalize the 
area under the theoretical curve, obtained when the full model of Fig. [5] is considered with our best fit parameters 
C5 (0) = 0.867, Maa = 0.985 GeV, to that under the experimental data. Here again it can be appreciated that, and 
despite the fact that the isospin factor of the A— pole mechanism (with the excitation of the A ++ resonance and its 
subsequent decay) in this channel is bigger than in the others, the effect of the background terms is quite significant. 

Next we show in Fig. \5\ the total ANL v^p — ■> /i~p7r + , v^n — > [i~pir° and v^n — ► n~nir + cross sections, with 
the kinematical cut W < 1.4 GeV, as a function of the neutrino energy and the predictions of the three schemes 
defined above: only AP contribution with C5 (0) = 1.2, Maa — 1-05 GeV and the full model derived in this work, 
including background terms, with the latter set of parameters for C^(q 2 ) and with that given in Eq. (f80|) . As can 
be appreciated in the different plots of the figure we achieve a reasonable description of data, finding the largest 
discrepancies in the ir + n channel. The inclusion of the chiral symmetry background terms derived in this work, brings 
in an overall improved description of the three channels as compared to the case where only the AP mechanism is 



16 



1.2 
1 



v^p -> fi pn + 



Full Model, Eq.(80) - 
Full Model, Eq. 48 

Only Delta, Eq.(48j ■ 

ANL I-K-H 
BNL (no 7rN cut) l — ■ — l 




0.3 
0.25 
^ 0.2 
0.15 
b 0.1 
0.05 




1 1 1 

Full Model, Eq.(80) 

■ Full Model Eq. 48 
Only Delta, Eq.(48) 

ANL i-H-H 
BNL (no 7rN cut) 1— ■— l 



0.35 
0.3 
0.25 



3 0.2 - 



"T 



Full Model, Eq.(80) - 

Full Model, Eq. 48 

Only Delta, Eq.(48) 

ANL I — ^ — I 

BNL (no 7rN cut) , a , H 




0.2 0.4 0.6 0.8 1 

E (GeV) 



0.6 0.8 1 

E (GeV) 



FIG. 5: Experimental and theoretical results for v^p — > fi~pn + , u^n — > fi~pir° and v^n — > fi~niT + cross sections, as a function of the 
neutrino energy. The ANL results |2(I and theoretical cross sections incorporate the kinematical cut W < 1.4 GeV. Dashed lines stand 
for the contribution of the excitation of the A resonance and its subsequent decay (AP mechanism) with C^(0) = 1.2 and Maa = 1-05 
GeV. Dashed-dotted and central solid lines are obtained when the full model of Fig. |2]is considered with cf(0) = 1.2, M AA = 1.05 GeV 
(dashed-dotted) and with our best fit parameters C^(0) = 0.867, Maa = 0.985 GeV (solid). In addition, we also show the 68% CL bands 
(solid lines) deduced from the Gaussian correlated errors quoted in Eq. H80I I. We also display BNL cross section data from Ref. :21] which 
do not include the W < 1.4 GeV cut in the irN invariant mass (see text). 



considered. In the case of Tr + n and ir°p cross sections the reduction of the contribution of this latter mechanism is 
compensated by the inclusion of the background terms. Our results are similar in quality to those obtained from the 
model of Ref. [l5|. We also display in the various plots of this figure, BNL cross section data from Ref. [21[ which 
do not include the W < 1.4 GeV cut in the irN invariant mass. For neutrino energies below 1 GeV, the effects of 
the wN invariant mass cut is almost negligible (see Table III of Ref. Such effects become much more sizeable 

for larger neutrino energies (see Table III of Ref. 20j) which have prevented us to present BNL cross section data 
in the plots for neutrino energies above 1 GeV. We observe some degree of inconsistency among the ANL and BNL 
measurements. The present model, including non-resonant background terms, with a C§{q 2 ) form factor consistent 
with the off diagonal Goldberger-Treiman relation (Eq. (|42p) would lead to a better description of the BNL data (see 
dashed-dotted lines). 

In Fig. [5] we compare the pion azimuthal and pion-nucleon invariant mass distributions (neutrino flux averaged) 
predicted by the different models examined here with that measured in the ANL experiment [20l. l26j| . For both plots, 
we normalize the area under the theoretical curve obtained when the full model of Fig.[5]is considered with parameters 
C^(0) = 0.867, Maa — 0.985 GeV, to that under the experimental data. The inclusion of chiral background terms 
leads to a more pronounced 0* dependence improving in this way the agreement with the observed event distribution 
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FIG. 6: Flux averaged ANL distribution of events in the pion azimuthal angle with W < 1.4 GeV (left) and in the ttN invariant mass 
(right) for /i _ p7r + final state. Data taken from Refs. |2(J and |2(I . respectively. Dashed lines stand for the contribution of the excitation 
of the A++ resonance and its subsequent decay (AP mechanism) with Cg (0) = 1.2 and Maa = 1.05 GeV. Dashed— dotted and central 
solid lines are obtained when the full model of Fig. [2] is considered with (7^(0) = 1.2, Maa = 1-05 GeV (dashed-dotted) and with our 
best fit parameters C^(0) = 0.867, Maa = 0.985 GeV (solid). In addition, we also show the 68% CL bands (solid lines) deduced from 
the Gaussian correlated errors quoted in Eq. 180> . 




FIG. 7: Energy dependence of the muon antincutrino cross section, with the invariant mass cut W < 1.4 GeV, for the P^n — > /j+7r _ n 
(left panel) and the tr~ exclusive production (right panel) reactions. Data— points are taken from the CERN PS experiment of Ref. |39| . 
Dashed lines stand for the contribution of the excitation of the A resonance and its subsequent decay (AP mechanism) with C^(0) = 1.2 
and Maa = 1-05 GeV. Dashed— dotted and central solid lines are obtained when the full model of Fig. [2] is considered with Cg^O) = 
1.2, M A A = 1-05 GeV (dashed-dotted) and with our best fit parameters C^{0) = 0.867, M A A = 0.985 GeV (solid). In addition, we also 
show the 68% CL bands (solid lines) deduced from the Gaussian correlated errors quoted in Eq. 180> . 



in the ANL experiment. In the right panel of Fig. [S] we show the W— distribution of ANL events, which clearly shows 
the A(1232) peak. The chiral background terms dominate the distribution near the pion production threshold, and 
they also produce a slight shift of the maximum of the distribution to lower invariant masses. 

In Fig. we compare the predictions of our model with the CERN-PS muon antineutrino cross section data of 
Ref. [39( | . Our model provides larger cross sections than the experiment, but nevertheless, we find here again a 
reasonable description of the data, which is certainly better than when the AP mechanism alone is considered. We 
would like to remind here that from isospin symmetry we have (rnr~ \j^ c _(0)\n) — (p7r + |j^ c+ (0)|p). Therefore the only 
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dynamical difference between the v^p — > /i p7r + , which we describe correctly (see Fig. [5]), and the v^n — > p, + pir 
reactions is the sign of the antisymmetric term of the lepton tensor (1$ = L^J. )■ Thus, if one neglects the parity- 
violating part of the antisymmetric hadronic tensor {Wqq^) ( see (ED m Appendix [A")) , the cross sections for 
both reactions will be the same but for differences in the interference between the axial and vector contributions of 
the hadronic current. However, this vector-axial interference does not affect the sum of the cross sections aiy^y — > 
fi~pTt + ) +a{D^n — > /i + p7r~), except for its contribution to the parity-violating part of the symmetric hadronic tensor. 
For instance at E — 3 GeV, the experimental value for this sum of cross sections is around 0.78 x 10 -38 cm 2 while 
our best theoretical prediction for it is around 11% bigger. Considering the large errors in experimental data we think 
this difference is not significative 14 . That might suggest that the discrepancy in the antineutrino cross section might 
be solved by the inclusion of relative phases between the vector and the axial current theoretical contributions. 

Relative signs among the different background contributions are well established, since all of them have been deduced 
from the same lagrangian (Eq. (|26|) ) and the vector and axial currents derived from it. One might think in a possible 
inconsistency between the relative signs of background and resonant pieces. However, 

• These relative signs are consistent with those deduced from a quark model picture [42| . 

• As already mentioned, the vector part of the model presented here, reduces to that derived in Ref. [3l[ for the 
eiV — * e'NTT reaction. Note the relation 

Cn<l 2 ) = ^JlMl 2 )/e (81) 

with e = y/4wa the proton charge, between CX and the usual 7 TV A coupling, / 7 , used in pion clectroproduction 
reactions (see Eqs.(5) and (A-12) in Ref. [31(). This relation is obtained from the non-relativistic reduction 
of the C\ Dirac structure, which dominates in that limit. Using f^jO) 0.122 as in Ref. [31| . one obtains 
C^(0) « 2.2 in good agreement with Eq. (|46|) . The model of Ref [3l| described 15 reasonably well the available 
data for the eN — > e'Nn reaction at low and intermediate energies, including pion angular dependences. This 
makes us confident on our election of relative signs between resonant and non-resonant terms. 

• In what the axial part concerns, we use a consistent sign convention for both the diagonal and off-diagonal 
Goldberger-Treiman's relations (see the choice of relative signs in Eqs. (|27|) and (|42l) ). 

Nevertheless, in the Appendix [B] we have examined the effect of including relative minus signs between the axial and 
vector resonant contributions and also between the AP and the background terms. Changing the relative sign between 
the vector and axial contributions of the A mechanism is totally discarded by the data, while modifying the relative 
sign between resonant and non-resonant terms has a little effect, once the W integration is performed. This can be 
understood by looking at the irN invariant mass distribution of the right panel in Fig. [6l At the A peak resonant 
and non-resonant contributions do not interfere, since the first one is purely imaginary while the latter one is real, 
and thus in this region the relative sign between both type of contributions is irrelevant. At lower and higher values 
of W, where the resonant contribution takes also a real part and thus it has a non-vanishing interference with the 
background terms, there exists a constructive and destructive, respectively, interference. A change of the relative sign 
between resonant and non-resonant terms would reverse constructive into destructive interferences, and vice- versa, 
but the net effect after integrating in W is greatly diminished. 

After this discussion, we stress that one would need relative phases (and not merely minus signs) between vector 
and axial resonant and non resonant contributions to improve the combined description of neutrino and antineutrino 
cross sections. However, very recently S.K. Singh and collaborators have pointed out [44| that nuclear medium effects 
and pion absorption (note that the antineutrino data of Ref. [39j presented in Fig. [7] were measured on freon-propane) 
were not properly taken into account in the original work of Bolognese and collaborators (39|. In this manner, the 
apparent discrepancies highlighted by the former discussion might disappear (see Figs. 17 and 18 of Ref. [3]). More 
accurate data, possibly available in the next future from the MiniBoonc and T2K experiments, in conjunction with 
Watson's theorem Q might shed light into this interesting issue. 

Next we pay attention to the differential cross section decomposition of Eq. (f74|) following the different allowed 
dependences on the pion azimuthal angle. In Figs. [5] and [5] we simultaneously compare results for the v^n — > ^pir a 
and v^p —> p, + nir° reactions. Thanks to isospin symmetry (see Eq. (|13|) ). the hadronic tensor is the same for both 



We note the discrepancy is significantly larger (rj 60%) if one uses C^(0) = 1.2. 

A small relative phase between the AP and the background terms was also included in Ref. [3l| . This phase is deduced from Watson's 
theorem 43], it depends on q 2 and W , and for the kinematics of interest in this work is comprised in the range 10 — 20°. 
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ANL Our results 



R+ = <r(yp — > vmx + ) / a(vp — » (j,~pir + ) 
Rq — a(vp — » i>pn°) / a(vp — ► fj,~pir + ) 
R~ — a(yn —> vpn~) / a(vp — > (i~pTr + ) 



0.12 ± 0.04 [49] 0.12 - 0.10 
0.09 ± 0.05 [49] 0.18 - 0.14 
0.11 ± 0.022 [48] 0.12 - 0.09 




TABLE I: NC to CC (up — > /i~piT + ) cross section ratios. Experimental data taken from the ANL analyses of Refs. [4g] (R-), 

and [4S| (Ro,R+). Our results are obtained for an incoming neutrino energy range of E = 0.6 — 1.2 GeV using our full model 
of Fig. [2] and with our best fit parameters C5 (0) = 0.867, Maa = 0.985 GeV. No kinematical cut in the W invariant mass has 
been used. 

processes and hence, they are only distinguished by the different neutrino- and antineutrino-induced lepton vertices, 
which produces a change of sign in the antisymmetric part of the leptonic tensor. Therefore, and following Eqs. (|A8[) 
and (|A9|) . the C* and E* structure functions are equal for both reactions, while the antisymmetric contributions to the 
neutrino and antineutrino A* , B* and D* ones change sign. We fix E = 1.5 GeV and q 2 — —0.5 GeV 2 , which naturally 
lies into the ANL kinematics, and we consider two different 7riV invariant masses, W = 1150 MeV and W = Ma, to 
better understand the effect of the chiral symmetry background terms on the structure functions defined in Eq. (|74[) . 
Results displayed in these plots show clearly both the difference between neutrino and antineutrino structure functions 
and the effect of the the chiral symmetry background terms on them. For instance in the W = Ma case, neutrino 
and antineutrino A* structure functions differ by about a factor of three, which will provide a similar factor in the 
integrated cross sections. In these channels, background terms have a greater influence in the antineutrino induced 
process than in the neutrino one. On the other hand, for W = 1150 MeV the AP mechanism becomes sub-dominant 
and the bulk of the structure functions is determined by the background terms and their interferences with the AP 
one. 

Besides, the interference between the AP mechanism and the rest of the chiral background terms leads to non- 
vanishing D* and E* structure functions. These functions provides dependences in sin0* and sin 20* and arise from 
the parity violating terms in the hadronic tensor decomposition in Eq. (|Aip . These parity- violating contributions to the 
fifth differential cross section d 5 <r Ul i / dfl(k')dE' 'c2f2(fc ff ) disappear when the pion solid angle integration is performed, 
as required by the scalar, up to the factor |, nature of the triple differential cross section d 3 'a Vl i / 'dVL(k')dE' '. 

Note that the coordinate system used to define dCl(k\) involves the pseudo-vector k x k' to set up the Y— axis, which 
induces the non-parity invariant nature of d 5 a/dfl(k')dE'dQ(k v ). In electropion production processes, the leptonic 
tensor is purely symmetric, and the symmetric part of the hadronic one can not contain terms involving the Levi-Civita 
tensor, since the electromagnetic interaction preserves parity. Hence in that case d 5 a / dQ(k')dE' dfl(k 7T ) turns out to 
be a scalar under parity. Moreover, these terms also induce T— odd correlations in the (L^ v,v >) illV W ,J ' v contraction (see 
Eqs. (|A1[) - (|A7|) in the Appendix [A"|) . which do not imply a g enuine violation of time reversal invariance because of 
the existence of strong final state interaction effects 16 [45|, |46[ . 



There hardly exist [47ll48ll49j] NC experimental measurements at intermediate energies. In the Gargamelle propane- 
freon experiment run at CERN [47|, the NC neutrino induced pion production cross sections, in all possible channels, 
were measured at averaged neutrino energy of around 2.2 GeV, and given in the form of NC/CC ratios. These 
data have been reanalyzed and absolute cross sections, without imposing any cut in the pion-nucleon invariant 
mass, have been recently provided [5fj| . Experiments using the Argonne 12-foot deuterium bubble chamber were run 
over a neutrino energy interval (0.3 < E < 1.5) GeV [H, |4^|. These latter experiments gave results for the NC 
vn — > vpir~ cross section [48|, as well as the NC to CC cross section ratios R+ = a(vv — > vnis + ) / <r(vp — > /i _ p7r + ), 
Po = a (yip — > h>pn a ) / 1 o(vp — > fi~pn + ) [49l ] and P_ = <y(yn — > vpn~) / ' a(vp — > /i~j>7r + ) [48]. 

The NC pion production reaction was proposed in Ref. [321 ] as a potential tool to distinguish r— neutrinos from 
antincutrinos, below the r— production threshold, but above the pion production one. Model independent neutrino- 
antineutrino asymmetries both in the totally integrated cross sections and in the pion azimuthal differential distri- 
butions were predicted in Ref. [32| . Results of the current model for these neutrino-antineutrino asymmetries were 
presented there. Here, we will focus on the comparison of our results with the available experimental data as well as 
emphasizing other aspects of the NC pion production processes. 



Within our formalism, the inclusion of the A resonance width accounts partially for the strong final state effects. 




B. NC pion production cross sections 
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[50] Our results 



a(up -f vpn°) 0.130 ± 0.020 0.105±0.006 

a(up -f vnw + ) 0.080 ± 0.020 0.091±0.003 

a(yn -» i/nvr ) 0.080 ± 0.020 0.104±0.006 

a(un -> up%-) 0.110 ± 0.030 0.082±0.003 

TABLE II: NC cross sections in units of 10 _38 cm 2 for different channels. Data shown correspond to the results of a recent 
reanalysis [5(j of the original data taken by the Gargamelle experiment at CERN [47j . Our results are obtained for an incoming 
neutrino energy of E — 2.2 GeV using our full model of Fig.[5]and with our best fit parameters C$(0) = 0.867 ±0.075, Maa = 
0.985 ± 0.082 GeV. No cut in the pion-nucleon invariant mass W has been applied and we Monte Carlo propagate the latter 
errors to our results for the cross sections. 



In Table [T] we compare our results for the P+ , Rq and i?_ NC over CC ratios with the ANL experimental data. 
Our results are obtained for an incoming neutrino energy range of E = 0.6 — 1.2 GeV using our full model of Fig. [3] 
and with our best fit parameters C^(0) = 0.867, Maa — 0.985 GeV. We find a fair agreement for all ratios, when the 
experimental uncertainties are taken into account. 

In Table |TT] we present results for the total NC cross sections in different channels at E =2.2 GeV and compare 
them with the reanalysis done in Ref. of the original data by the Gargamelle experiment at CERN ■ We use 
our full model of Fig. [2] with our best fit parameters C^(0) = 0.867, Maa = 0.985 GeV with no upper limit in the 
W— integration to allow a direct comparison with the experiment. The agreement with data is good. However a word 
of caution is in order here. As we discussed above, our model suffers from larger uncertainties in the W > 1.4 GeV 
region which will be accessible for neutrinos of this energy. 

In the left panel of Fig. [10] we show our results for the vn — > i/pw~ cross section as a function of the energy and 
compare them with the ANL data of Ref. [48[. There, and to better compare with the data, we also give results without 
the W < 1.4 GeV constraint. Up to incoming neutrino LAB energies of the order of 1 GeV the implementation of 
this cut hardly changes the cross section. We find a good description of the data. In the right panel of Fig. [10] we 
show the W— differential cross section for the same channel and neutrinos of 1 GeV. There one can appreciate clearly 
the A(1232) peak. The chiral background terms dominate the differential cross section near the pion production 
threshold, and they also produce a slight shift of the maximum of the cross section to lower invariant masses, as it 
also happened in the distribution of CC events displayed in Fig. [6] In both plots we see that the reduction of the 
contribution of the AP mechanism, due to the diminution of the value of C^(0), is partially compensated by the 
inclusion of the background terms. 

Next, we study the effect of the strange content of the nucleon within our model (left panel of Fig. fTTj) . We find that 
effects are even smaller than the statistical uncertainties deduced from the fit of the C§(q 2 ) form factor to the ANL 
data in Eq. (|50]) . Similarly, results of the right panel of this figure show that the isovector part of the NC completely 
dominates the pion production reaction at the intermediate energies studied here. This is because, as can be deduced 
from Eqs. (jf>4"]) . (|65[) and ([70]) . both cross sections just differ in the interference between the isovector and the isoscalar 
parts of the NC. Finally, we would like to mention that antineutrino induced cross sections are around a factor 2 or 3 
smaller than neutrino induced ones, as can be appreciated by comparing predictions for pn~ final state given in the 
left panel of Fig. [TU] and the right panel of Fig. [11] 



V. CONCLUSIONS 



We have developed a model for the weak pion production off the nucleon driven both by CC and NC at intermediate 
energies, which improves most of the existing ones. Besides the AP mechanism, we have also included some background 
terms, required by chiral symmetry. Starting from a SU(2) non-linear a model involving pions and nucleons, which 
implements the pattern of spontaneous chiral symmetry breaking of QCD, we derive the corresponding vector and 
axial currents (up to order 0(1//^))) which determine the structure of the chiral non-resonant terms. Vector current 
conservation and PCAC are also employed to establish some relations between the weak form factors. In this way 
constructed, this model represents the natural extension of that developed in Ref. [3l[ for the eN — > e'N'ir reaction. 

As a result of the inclusion of the background contributions, we had to re-fit the C$(q 2 ) form factor to the flux 
averaged v^p — > n~pn + ANL q 2 ~ differential cross section data with W < 1.4, finding a smaller contribution of the 
AP mechanism than traditionally assumed in the literature (see Eq. ([80]) ). We find a correction of the order of 
30% to the off diagonal Goldberger-Treiman relation (Eq. f4"2"]) L which we interpret is not in conflict with the lattice 
QCD results shown in Figure 4 of Ref. [30(, if they were extrapolated to realistic pion masses. Within this scheme, 
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we have calculated several differential and integrated 1 cross sections, including pion angular distributions, induced 
by neutrinos and antineutrinos and driven both by CC and NC. In all cases we find that the background terms 
produce quite significant effects, and for those quantities where there are experimental measurements, we find that 
the inclusion of these terms brings in an overall improved description of the data, as compared to the case where only 
the AP mechanism is considered. We give 68% CL bands for most of the computed observables as deduced from 
the Gaussian correlated errors quoted in Eq. (f80|) . For NC reactions the isoscalar contribution is quite small and in 
particular we find the nucleon strange content effects are smaller than the statistical uncertainties deduced from the 
fit of C^(q 2 ) to the ANL data. 

At higher irN invariant masses than those considered in this work, heavier resonances than the A(1232) (as for 
example iV(1440), AT(1535), iV(1520) • • • ) will certainly play an important role. However, we might safely expect 
the contribution of these heavier resonances to be negligible at pion threshold, where the chiral background terms 
computed in this work are dominant, and that it would remain quite small up to the ttN invariant masses around 
1.3-1.4 GeV considered here 

We also show that the interference between the AP and the background terms produces parity-violating contribu- 
tions to the fifth differential cross section d 5 a vi i/ dft(k')dE' dil(k 7T ) , which are intimately linked to T— odd correlations 
in the {L^ v,v ') ^ W^" contraction. However, these latter correlations do not imply a genuine violation of time reversal 
invariance because of the existence of strong final state interaction effects. 

The extension of this work to the study of the weak two-pion production off the nucleon near threshold is natural 
and will be presented elsewhere [5l[ . 
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APPENDIX A: DEPENDENCE OF THE NEUTRINO DIFFERENTIAL CROSS SECTION ON THE 

OUTGOING PION AZIMUTHAL ANGLE 

The hadronic tensor is completely determined by up to a total of 19 Lorentz scalar and real, structure functions [32| 
Wi(q 2 , p ■ q, p ■ k n , K ■ q), 

= (w& Ji:i< pc + (W£5JS' PV (ai) 

where (for simplicity from now on we drop the CC7r and [y] labels in the notation of the hadronic tensor) 

(W^f s C = WigT + W 2P »p v + W 3 q»q v + W 4 k%k» + W 5 (q»p» + q v p») + W^K + <fK) 

+ W 7 {jfK+^K) (A2) 

+ W 10 (fe^^pV + fc^ l Q/37 fc>V) (A3) 



(WnV - WiifeV -q v p») + W 1 2{q»K-q v k%) + W 13 [^K-p v K) (A4) 
{WnT = W 14 e^ al3 p a q p + W 15 e^ a0 Pa k^ + W 16 e^ al3 q a k^ + Wvr (<A^ 7 OV - <Z^ 7 OV) 

+ Wis (j^ 7 A£pV - P^OV) + Wi9 (fc^ Q/37 fc>V - fc^ 7 fc>V) (A5) 

The tensor (W^) = {W^ v ) s + i {W^) a when contracted with the leptonic one, L^v i provides a pseudo-scalar 
quantity, i.e., such contraction is not invariant under a parity transformation. Indeed, under a parity transformation 



There are some inconsistencies among the ANL and BNL totally integrated cross section data (Fig. [5)l and we would obtain a better 
description of the BNL data with a C^(q 2 ) form factor consistent with the off diagonal Goldberger-Treiman relation. 
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we have, 



42 -» ( L ^) H . (^) PV - - (^) PV (A6) 



PC PC PC ft') 

whereas the tensor (^) r " = (W^) s + i(W"% transforms as (L^f ' . This explains the origin of the 
adopted labels PC and PV, which stand for parity violating and conserving contributions to the fifth differential 
cross section d 5 cr^ l i/dft(k')dE'dil(k 7T ), respectively. The triple differential cross section d 3 cr^ l i/dil(k')dE' is a scalar, 
up to the factor |fc'|/|fc |. Thus all parity- violating contributions must disappear when performing the pion solid 
angle integration. A further remark concerns about the time-reversal (T) violation effects apparently encoded in the 
decomposition of the hadronic tensor in Eqs. (|A1HA5|) . Under a time reversal transformation, and taking into account 
the antiunitary character of the T— operator, we have 

£^(^) W , (^)pc^(n PC , (^r^-(n PV (at) 

and therefore (L^) is not T —invariant either, because of the presence of the PV terms in the hadronic tensor . 

This does not necessarily means that there exists a violation of T— invariance in the process because of the existence 
of strong final state interaction effects [45j, |46| . 

With our election of kinematics (k, k' in the XZ plane), we find (La )oy = (Ls)xy = [L^) Z y = (L^)ox = 

(4 )oa = (L^)xz = 0, and then 



L 



+ °° rlh f.2 



rlh h 2 f + oc jt, ul , 

^ (L^)„uWr = J ^[(L^UWT + 2(L^)o x W°* + 2(lM) te ie + (LW) 9X W? 
+ (4 v) )yyWr + {LW)„W? + 2(L^) xz Wr) 



= A s + B s cos 4> v + C s cos 2(f>„ + D s sin W 4- E s sin 2^ (A8) 

+°° rlh h 2 f + °° rlh h 2 r -i 

^pL (Li^wr = 2 / ^ (4'»)o»ie + (L^hyW^y + (L^) yz Wr\ 

= -A a - -BqCOS^ - D Q sin0 w (A9) 

where we explicitly show the </>„■ dependence 19 . The P V term of the hadronic tensor has led to the parity violating 
sin 0^ and sin 20^ contributions (all of them proportional to k vy ). They disappear when the pion solid angle integration 
is performed, as anticipated. The above equations automatically imply Eq. ([7]). 



APPENDIX B: VECTOR-AXIAL RESONANT AND RESONANT-NONRESONANT RELATIVE SIGNS 

Because of the poor description of the muon antineutrino cross section data of Ref. [3^] (see Fig. [Jj achieved with 
the different models studied in this work, we examined here the effect of including relative minus signs between the 
axial and vector resonant contributions and also between the AP and the background terms. We have focused on 
the flux averaged v^p — ► n~pn + ANL q 2 — differential cross section displayed in the left panel of Fig. QJ and on the 
total v^p — > ^i~pn + and v^n — > [i + pTt~ cross sections given as a function of the neutrino/antineutrino energy in 
the top and left panels of Figs. O of [JJ respectively. Results of our analysis are presented in Fig. [TJJ There, and 
in addition to the results obtained from the three models presented up to now (excitation of the A resonance and 
its subsequent decay, AP mechanism, with Cg (0) = 1.2 and Maa = 1-05 GeV (blue dashed-lines), and the full 
model of Fig. [2] with C£(0) = 1.2, Maa — 1-05 GeV (cyan dashed-dotted lines) and with our best fit parameters 
C*^(0) = 0.867, Maa = 0.985 GeV (red solid lines)) with the set of signs deduced in our scheme, we also show results 
from different choices of the relative signs. Curves denoted by V (A) have been obtained by changing the sign of the 
WNA vector (axial) form factors in Eq. (|40j) . and similarly results denoted by VA have been obtained by including 
an extra minus sign between the AP and the background contributions. Changing the sign of either the vector or the 
axial contributions of the A mechanism is strongly disfavored by the data, while modifying the relative sign between 
resonant and non-resonant terms has a little effect, as we already mentioned in the discussion of Fig. [7J 



Note that transformations given in Eqs. JB}, ||A6| | and iA7l imply that the leptonic tensor, by itself, is invariant under CPT. 

All structure functions, Wi— 1 ...19, depend on the Lorentz scalar p ■ k^ and k„ ■ q factors, which are functions of the angle formed between 

the q and k n vectors, and thus they are independent of </>„-, when q is taken along the Z— axis. 
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These results provide a further confirmation of the sign convention used in this work, and the muon antineutrino 
discrepancies of Fig. [7] point out to the existence of non-trivial relative phases (and not merely minus signs) between 
vector and axial resonant and non resonant contributions and/or more likely that some nuclear medium effects were 
not properly discounted [44j in Ref . (39j , when providing cross sections off the nucleon from measurements obtained 
in the liquid bubble chamber Gargamelle that was filled with propane and a small admixture of heavy freon CFsBr. 
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FIG. 8: Pion polar angle dependence of the 0* structure functions defined in Eq. H74H for the v^n — > n~pir° and u^p — > fi + nn° reactions. 
The neutrino incoming energy is E = 1.5 GeV, q 2 = —0.5 GeV 2 , W = Ma and the pion polar angle is referred to the nN CM frame. 
Neutrino (antineutrino) A*, B* and D* structure functions are displayed in the three upper right (left) hand side plots. The C* and E* 
structure functions, which are equal for both reactions, are showed in the two panels of the last row. Dashed lines stand for the contribution 
of the excitation of the A resonance and its subsequent decay (AP mechanism) with C^(0) = 1.2 and = 1.05 GeV. Dashed-dotted 

and central solid lines are obtained when the full model of Fig. [5] is considered with Cg (0) = 1.2, = 1-05 GeV (dashed-dotted) and 

with our best fit parameters C^-(O) = 0.867, Maa = 0.985 GeV (solid). In addition, we also show the 68% CL bands (solid lines) deduced 
from the Gaussian correlated errors quoted in Eq. 180> . 
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FIG. 10: Total (left) and W— differential (right) cross sections for the vn — > vpw~ reaction. Left panel data are taken from Ref. [48[ 
and have been measured without limiting the pion-nucleon invariant mass W . Dashed lines stand for the contribution of the excitation 
of the A resonance and its subsequent decay (AP mechanism) with C^(0) = 1.2 and = 1.05 GeV. Dashed— dotted and central solid 

lines are obtained when the full model of Fig. [5] is considered with C^(O) = 1.2, = 1.05 GeV (dashed-dotted) and with our best 

fit parameters C^(0) = 0.867, Maa = 0.985 GeV (solid). In addition, we also show the 68% CL bands (solid lines) deduced from the 
Gaussian correlated errors quoted in Eq. (1806 - In all these cases, we have limited the invariant mass phase-space (W < 1.4 GeV). In the 
left hand side panel the long dashed-line stands for our full model results without including the W < 1.4 GeV cut. The LAB incoming 
neutrino energy in the right hand side plot is 1 GeV. 



o.os 




E (GeV) 



0.04 
0.035 

0.03 
0.025 

0.02 
0.015 

0.01 
0.005 






1 


1 1 






vn - 








vp - 


■* pniT + 






^ i 


1 1 





0.5 



1.5 



E (GeV) 



FIG. 11: Total up — > upir (left panel), vn — * vpir~ and Dp — » Pn7r+ (right panel) cross sections, with the W < 1.4 GeV cut, as a function 
of the neutrino or antineutrino energy. In the left panel we show results obtained for two different values of the nucleon strange content 
(gs in Eq. I I72H ) with the full model of Fig.[2]and our fitted C^(q 2 ) form factor. We also include the 68% CL bands inferred from Eq. (180 6 . 
In the right panel we do not give the CL bands and gs is set to —0.15. 



27 



1.2 



o 
i — i 

b 



lb 




1 1 

Delta 

Full Eq.(48) 

Full Eq. 80) 

Delta. V 

Full Eq.(48) V 

Full Eq. 80 V 

Full Eq. 48 A 

Full Eq. 80} A 

Full Eq.(48) VA 

Full Eq.(80 VA 

ANL H-K- 



Delta 

Full Eq.(80) 

Full Eq. 48} V 

Full Eq. 48} A 

Full Eq.(48) VA 

l — X — l ANL 




Full Eq.(48) 

Delta V, 

Full Eq.(80 V 

Full Eq. 80 A 

Full Eq.{80 VA 

i-K-h CERN-PS 
n 1 1 1 r 




E (GeV) 



E (GeV) 



FIG. 12: Flux averaged v^p — > fj,~pir + ANL q 2 — differential cross section (top) and total u^p — > fj,~pir + and P M n — > )i Jr p-K~ cross sections 
as a function of the neutrino/antineutrino energy (bottom) from different choices of the vector— axial resonant and resonant— nonresonant 
relative signs (see text for details). 



